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ABSTRACT 

A  complete  set  of  thermal  conditions,  mechanical  loads,  geometry  as 
well  as  thermal,  elastic,  fast  fracture,  and  time  dependent  material  pro¬ 
perties  was  furnished  by  ANMRC  for  the  hub  of  a  gas  turbine  disc.  A  finite 
element  computer  model  was  prepared  for  the  disc  from  this  data.  The  tem- 
perture  and  stress  distributions,  the  fast  fracture  and  the  time  dependent 
reliabilities  were  calculated.  The  results  were  presented  in  iso-stress 
and  iso-thermal  plots  and  in  graphical  data. 
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Foreword 


This  report  represents  the  analysis  of  a  structural  ceramic  com¬ 
ponent  selected  by  AMMRC.  The  component  selected  was  the  hub  of  a  hot 
pressed  silicon  .  nitride  turbine  rotor.  The  required  geometrical, 
material,  strength,  and  time  dependent  data  was  supplied  by  AMMRC. 
This  work  was  conducted  as  part  of  the  "Methodology  for  Ceramic  Life 
Program,"  initiated  by  Dr.  Robert  Schulz  of  the  Office  of 
Conservation,  division  of  Transportation  Systems,  Department  of 
Energy,  and  monitored  by  the  Army  Materials  and  Mechanics  Research 
Center  under  Contract  Number  DAAG-46-77-C-0028.  This  work  was  part  of 
the  continuing  investigation  of  analytical  and  experimental  methods  in 
ceramic  life  prediction  aimed  toward  utilizing  structural  ceramics  in 
high  temperature  applications.  The  principal  Investigator  of  this 
program  was  R.  R.  Baker,  Ceramic  Materials  Department,  engineering  and 
Research  Staff,  Ford  Motor  Company.  The  technical  monitor  was  Dr.  E. 

M.  Lenoe  of  AMMRC.  The  author  wishes  to  thank  Drs.  E.  M.  Lenoe  and  R. 

N.  Katz  of  AMMRC  for  suggestions  in  carrying  out  the  program. 
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1 .  INTRODUCTION 


This  report  presents  the  analysis  of  a  gas  turbine  disc  based  on 
data  supplied  by  AMMRC.  The  data  supplied  was  the  disc's  geometry, 
its  rotational  speed,  forces  to  simulate  the  blade  loads,  heat 
transfer  coefficients  and  corresponding  fluid  temperature.  Also 
supplied  were  the  material's  elastic,  fast  fracture,  and  time  dependent 
properties. 

The  primary  purpose  of  the  analysis  was  to  calculate  reliability 
versus  time  for  the  steady  state  speed  and  temperature  conditions. 
The  time  dependent  theory  used  is  explained  in  Section  2.  In  order  to 
complete  a  time  dependent  reliability  analysis  it  is  necessary  to 
calculate  the  reliability  at  time  equals  zero  which  is  the  fast  frac¬ 
ture  reliability.  Weibull  theory  was  used  for  calculation  of  fast 
fracture  reliability  and  the  theory  is  briefly  reviewed  in  Section  3. 

The  stresses  and  temperatures  were  calculated  with  finite  element 
computer  programs  and  the  results  are  plotted  in  contour  plots.  For 
reference  the  stresses  due  to  mechanical  loads  at  room  temperature 
were  calculated  and  plotted.  The  data  supplied  by  AMMRC  was  repro¬ 
duced  here  in  the  tabular  form  required  by  the  computer  codes  used  In 
the  analysis. 

2.  TIME  DEPENDENT  THEORY 

The  prediction  of  time  to  failure  in  structural  ceramics  is 
usually  based  on  the  equation 

V  -  AKTn  .  (1) 

This  equation  relates  the  velocity  of  a  crack  in  the  material  to  the 
stress  intensity  factor  of  the  tip  of  the  crack.  Sivems(l)  observed 
this  behavior  in  steel  and  published  velocity  versus  stress  intensity 
data  and  fitted  the  constant  "A"  and  crack  velocity  exponent  "nH  to 
the  data  with  linear  regression  analysls(2).  EvansC^)  observed  simi¬ 
lar  behavior  in  ceramic  materials  and  using  the  fracture  mechanics 
equation 

Kt  -  oYa1/2  (2) 

arrived  at  the  following  equation  for  calculating  crack  growth  versus 
time(^) . 


Equation  3  applies  only  to  a  single  crack  in  a  uniform  stress  field. 
To  allow  calculation  of  time  dependent  failure  in  complex  structures. 
Paluszny(5)  combined  the  ideas  in  Equation  3  with  the  Weibull W 
equation  for  calculating  fast  fracture  reliability  in  brittle 
materials  and  arrived  at  Eqn.  4. 


R  ■  exp 
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Equation  4  was  reduced  to  Eqn.  5  In  Reference  7. 


n,. 


In  this  study  Eqn.  5  was  used  to  calculate  the  reliabilities  versus 
time  of  complex  structures  by  dividing  the  structure  Into  small  ele¬ 
ments  where  the  stress  was  considered  constant  and  evaluating  Eqn.  5, 
the  reliability  of  a  structure  with  HqM  elements  was  then  calculated 
with  Eqn.  6. 

R  m  R,  .  R-  .  R,  ...  R  (6) 

12  1  q 

This  technique  Is  suitable  for  use  with  finite  element  stress  programs 
where  Rj  Is  evaluated  for  each  element.  The  constant  In  Eqn.  5  Is 
given  as  Eqn.  7(5). 


(n-2)YZAKIcn~Z 

Equations  5,  6,  and  7  were  used  In  this  report  to  calculate  time 
dependent  reliability. 

3.  WE  I BULL  THEORY 

The  evaluation  of  Eqn.  5  requires  that  the  fast  fracture  reliabi¬ 
lity  be  evaluated.  For  this  study,  the  Welbull(&)  equation  was  used. 
Welbull  proposed  that  each  element  of  a  structure  has  a  definite  pro¬ 
bability  of  failure  and  that  the  entire  structure  could  be  considered 
to  be  an  assembly  of  the  individual  elements  and  their  associated  pro¬ 
bability  of  failure.  Vardar  and  Finnie(8)  give  an  integral  for¬ 
mulation  of  the  Welbull  approach: 


«  1  -  exp 


VK  I, 


dA)dV 

n 


The  term  in  the  parentheses  Is  evaluated  on  the  surface  of  the 
unit  sphere  (Fig.  1),  over  the  regions  where  the  normal  stress  Is  ten¬ 
sile  and  neglecting  regions  where  the  normal  stress  is  compressive. 
The  reasons  for  neglecting  compressive  stress  are  1)  for  structural 
ceramics,  compression  Is  not  as  detrimental  In  causing  fracture  as  Is 
tension,  and  2)  mathematical  complexities  are  minlmlzed(9) .  The  nor¬ 
mal  stress  on  the  surface  of  the  unit  sphere  In  terms  of  the  maximum 
principal  stresses  and  the  polar  and  azimuthal  angles  is  given  by: 


2. 


dA=cos  ^d^d  ^ 


Fig.  1.  Geometric  variables  describing  location  on  unit  sphere 

2  2  2  2 
®n  «  cos  ^(OjCos  »|i  +  o2sin  \|i)  +  OjSin  <fr  (9) 

The  constant,  K,  is  given  by: 


K 


2m  +  1 

_  m 

2iro 

o 


(10) 


The  (2m+l)/2ir  term  in  Eq.  10  is  a  compatability  factor  required  to 
make  the  result  of  integrating  Eq.  8  using  the  normal  stress  distribu¬ 
tion  of  Eq.  (9)  for  uniaxial  stress  cases,  agree  with  the  results 
obtained  from  the  one-dimensional  Weibull  equation: 


*f  -  *  "  «P  dV  (11) 
o 

The  integration  of  Eq.  8  over  an  entire  structure  is  accomplished 
by  dividing  the  structure  into  finite  elements  in  which  the  principal 
stresses  can  be  assumed  to  be  constant.  The  result  of  integrating 
over  the  unit  sphere  is  constant  in  each  finite  element;  therefore, 
the  risk  of  rupture  for  each  element  is  simply  the  volume  of  the  ele¬ 
ment  times  the  Integrated  result  over  the  unit  sphere.  This  procedure 
is.  very  convenient  when  used  with  finite-element  stress  programs  and 
is  the  one  used  to  calculate  the  probabillty-of-failure  distributions 
for  the  structures  in  the  present  paper. 


li 
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In  the  Weibull  formulation,  o0  and  n  are  parameters  which  show 
the  probability  of  having  a  certain  strength  in  any  finite  element. 
The  values  of  oQ  and  m  are  obtained  from  experimental  test  results. 
Weibull  showed  that  the  characteristic  strength,  0o,  is  identical  to 
the  ultimate  strength  of  the  material  in  the  classical  theory  as  the 
Weibull  modulus,  m,  Increases  indefinitely.  This  result  provides  a 
physical  interpretation  of  the  characteristic  strength.  The  Weibull 
modulus  may  be  interpreted  as  a  measure  of  variability  in  the  strength 
distribution;  as  m  Increases,  the  variability  of  the  strength 
diminishes. 


For  this  report  value  of  o0  was  calculated  with  Eqn.  12(10). 


a 

o 


MORq  (-j£) 


1/m 


1/m 


(12) 


The  parameters  oD  and  m  define  the  strength  of  the  material  and 
are  Independent  of  the  volume  and  geometry  of  the  structure.  They  are 
the  values  used  to  determine  the  failure  distribution  of  a  complex 
structure  and  are  generally  determined  from  MOR  bar  tests.  It  is 
important  to  note  that  the  MOR  is  not  an  intrinsic  material  parameter 
but  depends  on  the  geometry  of  the  MOR  bar  and  test  fixture.  A  con¬ 
venient  formula  for  scaling  the  MOR  value  from  one  geometry  and  test 
fixture  to  another  is: 


morb 

MOR. 

A 


/LM+  mL*2> 


1/m 


(13) 


The  subscripts  here  identify  bar  and  fixture  A  and  bar  and  fix¬ 
ture  B.  This  formula  assumes  four-point  loading  with  the  inner  span 
place  symmetrically  with  respect  to  the  outer  span.  This  formula  was 
used  to  scale  the  MOR  values  furnished  by  AMMRC  to  standard  "A-size" 
MOR  values.  The  Weibull  equation  has  been  shown  to  be  very  effective 
for  predicting  the  experimental  failure  distributions  of  complex 
structure  from  data  bases  obtained  on  MOR  bar  tests(8),  (10). 


4.  FINITE  ELEMENT  MODEL 


The  finite  element  model  set  up  to  simulate  the  hub  is  shown  in 
Fig.  2.  It  has  216  elements  and  243  nodes.  The  elements  and  nodes 
are  numbered  horizontally  from  right  to  left,  starting  at  the  bore. 
Element  numbers  are  located  in  the  center  of  the  element,  and  node 
numbers  are  shown  on  the  exterior  of  the  grid.  Node  numbers  236,  237, 
238,  240,  241,  and  242  were  located  with  respect  to  the  centerline  to 
coincide  with  the  location  of  the  forces  used  to  simulate  the  blade 
loads.  Not  all  node  and  element  numbers  are  shown  on  the  mesh  in  Fig. 
2.,  but  any  one  may  be  located  by  counting  horizontally.  The  node  and 
element  numbers  are  needed  for  Interpreting  the  forces  used  to  simu¬ 
late  blade  loads  and  the  heat  transfer  data. 
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Fig.  2.  Disk  finite  element  grid. 


5.  MATERIAL  DATA 

The  calculation  of  temperatures,  stresses,  fast  fracture  and  time 
dependent  reliabilities  requires  a  number  of  material  properties. 
These  properties  were  supplied  by  AMMRC  or  it  was  suggested  that  a 
particular  property  was  similar  to  the  property  of  a  known  material. 
The  properties  used  are  given  in  the  Appendix  B.  AMMRC  suggested  that 
the  thermal  properties  of  Norton's  NC-132  be  used  and  thus  the  thermal 
properties  given  in  Table  I  are  for  NC-132. 

The  elastic  properties  are  given  in  Table  II.  Young's  modulus, 
Poisson's  ratio,  and  shear  modulus  were  furnished  in  graphical  form 
and  are  reproduced  here  in  tabular  form.  The  coefficient  of  thermal 
expansion  was  furnished  in  tabular  form,  Table  II. 

AMMRC  furnished  modulus  of  rupture  data,  giving  the  values  for 
50Z  probability  of  failure;  however,  the  size  of  the  MOR  bar  was  not 
that  of  a  standard  "A”  size  bar.  This  required  two  conversions  to  get 
the  MOR  bar  data  into  the  values  accepted  by  Ford's  codes  for  calcu¬ 
lating  fast  fracture  reliability.  First  the  characteristic  MOR  was 
calculated  from  Eqn.  14 


Pj  ■  1  -  exp  - 


(14) 


(MOR  )® 
MOR 


o 


where  Pf  -  0.5,  HOR  =  the  value  furnished  for  50%  probability  of 
failure,  and  m=15.  The  results  were  then  scaled  to  the  standard  "A" 
size  bar  using  Eqn.  13,  which  with  the  bar  dimensions  and  the  given 
Weibull  modulus  of  15  becomes 


M°Rb 

MOR, 

A 


1.0367 


(15) 


The  results  after  these  two  conversions  are  presented  in  Table  III. 

The  time  dependent  material  properties  are  given  in  Table  IV. 
The  values  were  furnished  for  1100,  1150,  1200,  and  1250°C.  Most  of 
the  computed  tempertures  in  the  disc  were  below  1100°C,  with  tem¬ 
peratures  down  to  787°C.  Since  values  of  B  and  n  are  required  at  all 
temperatures  the  B  and  n  values  furnished  for  1100°C  were  assumed  to 
apply  down  to  room  temperature.  The  alternate  approach  in  the  absence 
of  values  for  B  and  n  is  to  use  a  lower  limit  on  temperature  below 
which  no  crack  growth  Is  assumed  to  occur.  This  approach  was  not  used 
in  this  report,  since  a  lower  limit  was  not  furnished.  The  value  of  B 
furnished  was  assumed  to  be  the  one  obtained  from  the  application  of 
Eqn.  7. 


6.  MECHANICAL  STRESSES 


The  primary  purpose  of  this  report  is  to  calculate  time  dependent 
reliabilities  and  the  associated  temperatures  and  stresses;  however, 
for  interpreting  the  results  it  is  useful  to  have  the  stresses  due  to 
centrifugal  effects.  This  section  presents  these  stresses  at  room 
temperature.  Table  V  shows  the  forces  used  to  simulate  blade  loads 
and  the  node  number  at  which  the  particular  load  was  applied.  This 
data  was  furnished  in  tabular  form  by  AMMRC.  The  mesh  was  constructed 
so  that  the  nodes  were  locted  to  correspond  with  the  geometric  loca¬ 
tion.  The  forces  were  applied  radially  outward  at  these  nodes. 

The  rotational  speed  at  which  the  stresses  were  calculated  was 
61,500  rpm.  The  maximum  principal  stresses  ar-'  plotted  in  Fig.  3. 
The  bore  stress  is  165  MPa,  and  the  neck  stress  is  125  MPa.  The 
tangential  and  radial  stresses  are  plotted  in  Figs.  4  and  5  respec¬ 
tively.  The  isostress  plots  are  in  increments  of  5  MPa.  The  fast 
fracture  probability  of  failure  for  these  stresses  and  the  Weibull 
data  given  in  Table  III  is  3.725xl0”8. 

7.  COMBINED  THERMAL  AND  MECHANICAL  STRESSES 


This  section  presents  the  results  of  the  analysis  of  the  combined 
thermal  and  mechanical  stresses.  In  the  analysis  a  steady  state  tem¬ 
perature  distribution  and  no  creep  was  assumed.  The  heat  transfer 
data  furnished  was  in  graphical  form.  In  order  to  input  the  data  into 
the  computer  code  it  was  necessary  to  convert  the  data  into  numerical 
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form.  Table  VI  gives  the  data  as  prepared  for  the  computer.  Shown 
are  the  heat  transfer  coefficients  and  adiabatic  wall  temperatures  and 
the  corresponding  element  and  node  numbers.  Data  was  supplied  for  the 
pedestal,  and  upstream  and  downstream  faces  of  the  disc.  The  data 
furnished  assumed  no  heat  transfer  in  bore  and  at  the  attachments; 
therefore,  all  surfaces  not  listed  in  Table  VI  were  treated  as  being 
adiabatic.  The  steady  state  temperatures  calculated  from  these 
boundry  conditions  are  plotted  in  Fig.  6.  Temperatures  shown  are  in 
degrees  centigrade,  and  the  maximum  calculated  nodal  temperature  is 
1 121°C  at  node  216. 

The  combined  stresses  due  to  the  thermal  gradients  shown  in  Fig. 
6  and  the  centrifugal  loads  given  in  the  previous  section  were  calcu¬ 
lated  and  were  plotted  in  Figs.  7,  8,  and  9.  These  are  the  maximum 
principal,  tangential,  and  radial  stresses,  respectively.  The  maximum 
bore  stress  is  220  MPa  and  the  maximum  stress  in  the  neck  is  200  MPa. 
The  gradient  between  the  isostress  lines  is  5  MPa. 


Fig.  5.  Radial  stresses  (MPa]  in  disc,  centrifugal  stresses 


8. 


Fig.  6.  Isotherms  [°C]  in  disc  of  61500  rpm. 


Fig.  7.  Maximum  principal  stresses  {MPa]  in  disc,  combining  thermal 
and  centrifugal  stresses  at  61500  rpm. 
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8.  RELIABILITIES 


The  reliability  versus  time  for  the  disc  was  calculated  for  up  to 
1000  hours  of  steady  state  running  and  the  result  is  plotted  in  Fig* 
10.  The  fast  fracture  reliability,  that  is  the  reliability  at  time 
equal  zero,  is  0.99982.  The  time  dependent  parameters  used  are  in 
Table  IV.  The  results  assume  continuous  opertion  at  61500  rpm  at  the 
temperatures  shown  in  Fig.  6.  The  effects  of  creep  and  oxidation  are 
assumed  to  be  negligible.  The  maximum  principal  stress  in  the  bore 
225  MPa  at  a  temperature  of  823°C.  The  reliability  for  steady  state 
running  at  these  conditions  is  at  time  equal  zero,  0.99982,  and  at 
time  equal  1000  hours,  0.667. 


TIME  (HOURS) 

Fig.  10.  Reliability  versus  time  for  steady  state  operation. 

9.  SUMMARY 

The  mechanical  stresses  and  combined  mechanical  and  thermal 
stresses  were  calculated  for  a  hot  pressed  silicon  nitride  turbine 
disc.  At  room  temperature  the  maximum  principal  stress  was  165  MPa 
due  to  mechanical  loads  at  61500  rpm.  The  corresponding  fast  fracture 
probability  of  failure  was  3.725xl0~8.  The  steady  state  temperature 
distribution  was  calculated  for  the  boundary  conditions  at  61500  rpm. 
The  maximum  temperature  in  the  bore  was  824°C,  the  maximum  temperature 
in  the  neck  was  946°C,  and  the  maximum  temperature  on  the  platform  was 
1121°C.  The  thermal  stresses  due  to  this  temperature  distribution 
were  combined  with  the  stresses  due  to  mechanical  loads  at  elevated 
temperatures.  At  these  conditions  the  maximum  principal  stress  in  the 
bore  was  225  MPa  at  823°C,  and  the  maximum  principal  stress  in  the 
neck  was  210  MPa  at  929°C. 
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The  reliabilities  versus  time  were  calculated,  assuming  no  creep 
and  continuous  running  at  61500  rpm,  for  the  combined  stress  and 
steady  state  temperature  conditions.  At  time  equal  to  zero  the  calcu¬ 
lated  reliability  is  0.99982  and  at  1000  hours  the  reliability  is 
0.667. 
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11.  APPENDIX  A  NOMENCLATURE 

a  Crack  length,  m 

aQ  Initial  crack  length,  m 

2 

A  Area  on  unit  sphere,  m 

A  Premultlpller  In  crack  velocity  equation,  meter/second  (MPa-m)n 

2 

B  Constant  In  reliability  versus  time  equation,  MPa  -hr 

B  Risk  of  rupture 

b  Width  of  MOR  specimen,  m 
h  Height  of  MOR  specimen,  m 

K  Constant,  defined  by  Eqn.  10 

Kj  Stress  Intensity  MPav'm 

KfC  Critical  stress  Intensity,  a  material  parameter,  MPa/m 

Li  Bottom  span  of  MOR  specimen,  m 

L2  Top  span  of  MOR  specimen,  m 

m  Weibull  Modulus 

MOR  Modulus  of  rupture,  MPa 

MORq  Characteristic  modulus  of  rupture,  MPa 

n  Crack  velocity  exponent 

Pf  Probability  of  failure 

R  Reliability  as  a  function  of  time 

Rff  Fast  fracture  reliability 

o  Stress,  MPa 

og  Characteristic  strength  of  the  structure,  MPa 

On  Normal  stress,  MPa 

0l>°2t°3  Maximum  principal  stresses,  MPa 
oQ  Weibull  parameter,  MPa/(m3)l/® 
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t 


Time,  hours 

3 

V  Volume,  m 

V  Crack  velocity,  m/sec 

V  Stress  Intensity  factor  coefficient,  non-dimensional 

4>  Azimuthal  angle 

$  Polar  angle 


TABLE  I 


Hot  Pressed 
Nitride 
Density  ■  3 


Thermal  Properties 


Temperature  Thermal  Specific 

Conductivity  Heat 
°C  watt  joule 

_ meter  *C~ _ kg  *C 


Silicon 

21. 

29.4 

754 

260. 

26.0 

963 

19  gm/cm3 

538. 

22.5 

1088 

816. 

19.0 

1214 

TABLE  II 


Elastic  Properties 


Temp. 

Young ' s 
Modulus 

Poisson's 

Ratio 

Shear 

Modulus 

Hot  Pressed  Silicon 
Nitride 

°C 

GN/m2 

GN/m2 

Density 

23. 

320. 

.280 

122. 

■3.19  gm/cm3 

600. 

311. 

.278 

121. 

Coefficient 
of  Thermal 
Expansion 

1/°C 


1.13 

3.39 


TABLE  III 

Fast  Fracture  Properties 

Temperature 
°C 


Characteristic  Weibull 

MOR  Modulus 

MPa 


TABLE  IV 


Time  Dependent  Properties 

Temperature  Parameter  Crack  Velocity 

B  Exponent 


°c 

MPa2-Sec 

21 

515. 

40 

1100 

515. 

40 

1150 

515. 

40 

1200 

8.81 

40 

1250 

.12 

40 

19 


TABLE  V 


Forces  Simulating  Blade  Loads 

Node  Force 

Number  Newtons 


236 

5148.8 

237 

12821.9 

238 

15903.6 

240 

15996.9 

241 

13617.0 

242 

5631.7 

TABLE  VI 


Pedestal 


Upstream 
Face  of  Disc 


Downstream 
Face  of  Disc 


Heat  Transfer  Data 


Element  Node  Node 

Number  1  2 


Heat  Transfer 
Coefficient 
W 


Adiabatic 
Wall  Temperature 

°C 


216 

243 

242 

215 

242 

241 

214 

241 

240 

213 

240 

239 

212 

239 

238 

211 

238 

237 

210 

237 

236 

209 

236 

235 

100 

121 

120 

99 

120 

119 

110 

119 

130 

118 

130 

137 

124 

137 

144 

130 

144 

151 

136 

151 

158 

142 

158 

165 

148 

165 

172 

154 

17  2 

179 

160 

179 

186 

166 

186 

193 

172 

193 

200 

178 

200  ' 

207 

192 

207 

216 

200 

216 

225 

208 

225 

234 

216 

234 

243 

88 

95 

109 

101 

109 

122 

111 

122 

131 

119 

131 

138 

125 

138 

145 

145 

152 

137 

152 

159 

143 

159 

166 

149 

166 

173 

155 

173 

180 

161 

180 

187 

167 

187 

194 

173 

194 

201 

179 

201 

208 

193 

208 

217 

201 

217 

226 

209 

226 

235 

2650. 

1480. 

1260. 

1230. 

1200. 

1180. 

1150. 

1130. 

1170. 

1150. 

1090. 

1030. 

970. 

925. 

885. 

850. 

815. 

790. 

760. 

730. 

710. 

690. 

660. 

640. 

625. 

608. 

580. 

550. 

520. 

490. 

465. 

450. 

435. 

420. 

410. 

400. 

390. 

383. 

380. 

370. 

365. 

355. 

350. 


1250 

1250 

1250 

1250 

1250 

1250 

1250 

1250 

700 

701 
709 
712 
722 
730 
737 
745 
750 
757 
763 
772 
776 
795 
820 
900 

1042 

1150 

755 

757 

763 

776 

786 

795 

800 

810 

815 

825 

835 

845 

900 

995 

1100 

1231 
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